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We study the directional dependence of the angular two-point correlation 
function in maps of the cosmic microwave background (CMB). We propose two 
new statistics, one which measures the correlation of each point in the sky with 
a ring of points separated angle 9 away, and a second that measures the missing 
angular correlation above 60 degrees as a function of direction. Using these 
statistics, we find that most of the low power in cut-sky maps measured by the 
WMAP experiment comes from unusually low contributions from the directions 
of the lobes of the quadrupole and the octupole. These findings may aid a future 
explanation of why the CMB exhibits low power at large angular scales. 



1. Introduction 

With the launch of the Wilkinson Microwave Anisot ropy Probe (WMAP) in 2003, the 



CMB has been measured in highly detailed ful 



over t he past few years (ICruz et al. 



2005; Hoftuft et al. 2009 



20071 ; iHinshaw et al.l 120091; iJarosik et al.l 1201 11). which have been examined in great detail 



12005 



Hansen et al 



sky maps ([Bennett et al.ll2003l ; ISpergel et ai- 



de Oliveira-Costa et al 



2009 



Kim fc Naselsky 



2004 



2010 



Land &; Magueiio 



Bennett et al 



20111 ). 



In particular, the angular two-point correlation function C (6) has been studied; it is defined 
as the average product between the temperature of two points angle 6 apart 



C(9)=T(Q 1 ),T(Q, 



If2l-n 2 =cos(0)- 



Here T(Cl) is the fluctuation around the mean of the temperature in direction Cl on the 
sky. Several anomalies have been claimed in the angular co rrelation function , especially 
the missing power on large angular scales (for a review, see ICopi et al.l ( 120101 ) ) . Specifi- 
cally, the angular correlation function is very nearly zero at scales above 60°; such a low 
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correlation has a significance of > 3.2a in the standard Gaussian random, statistically 
isotropic cosmology. This re sult, first observed in COBE data, has been strengthened 
with WMAP first-year data (IHinshaw et al.lll996l ). as well as later WMAP observations 



JCopi et all [2007[ hood: Isarkar et al 



significance ( Pontzen fc Peirisl 201ol 



20111). though questions have been ra ised regarding its 



Efstathiou et al.ll20ld Ma et al.ll2011[ ) 



In this brief report, we explore the directional contributions to the angular two-point 
correlation function. Directional information is lost in the original definition of C(9), as 
Eq. (pQ) assumes statistical isotropy and averages over all directions on the sky. Our goal 
is to provide a generalization of Eq. (JTJ, and ascertain if there are specific directions which 
result in the unusually low correlation. 



2. Fixed-vertex C(9) 

The angular two-point correlation function is defined by Eq. (pQ), where the average is 
taken over all pairs of points separated by angle 9 on the sky. 

We now propose a new statistical quantity, a fixed-vertex correlation function 

C(Cl,9) = T(n),T(Cl 2 )\^ 2=cos(e) (2) 

where Cl is a direction in the sky, and the average is taken over a ring of pixels separated by 
9 from the direction Cl. In particular, notice that C(9) = jj J2i=i 9)!i Maps of C(0, 9) 
for 9 = 60° and 9 = 120° are shown in Figure 1. 

In order to quantify the low p ower in the angular two-point function between approxi- 
mately 60° and 180°, the statistic JSpergel et al.l l2Q03h 

f l/2 

Si/2= / C(9) 2 d(cos0) (4) 



1/2 

has been widely used. In order to examine possible anisotropies in C(9), we attempt to 



1 On a masked sky, the standard C(0) is no longer the simple average of C(Ct, 9). We need to introduce 
the weight function W(Q,,6), the fraction of unmasked points in the ring 9 degrees away from Cl. Thus, 



C(9) = 



(3) 



Qi ■Q,2—cos(0) 



Fig. 1. — Full-sky maps of the fixed- vertex correlation function C(Cl,9) for 9 = 60° (left 
panel) and 9 = 120° (right panel) at resolution N s ^c — 64, calculated for the cut-sky ILC 
map. 




Fig. 2. — The relationship between Sy 2 and S±/2 is shown. The first black curve (labeled 
'Total S1/2') indicates the median relationship in random Gaussian skies. The blue and pur- 
ple bands surrounding it are (respectively) 68% and 95% contours, indicating the spread of 
the distribution due to cosmic variance. (Note that the contours shown are only approximate 
for very high or low S1/2 values, due to very low statistics.) The second black curve (labeled) 
logarithmically plots the relative fraction of random maps; the green and red asterisks indi- 
cate the values for S®2 and Si/ 2 from the fiducial cut-sky and full-sky ILC map, respectively. 
The narrow contours indicate a very tight relation between the two statistics. The cut sky 
(green) fits the typical relation, but has an unusually low value for ; hence, it is possible 
that S1/2 is low because S®2 is. See text for details. 
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separate 61/2 into contributions from specific directions Cl; we define 

S 1/2 (£l) = 1 ' C(Cl,9)C(9)d(cos9), (5) 

so that the sum of 6*1/2 0) over all pixels is 6i/ 2 o 

Note that the quantity 6i/ 2 (f2) is non-local, in that it receives contributions from not 
only the direction Cl, but from the whole sky, due to the term C(9) in Eq. (JSJ). So, for 
example, excluding a pixel Q will not change the value of 5*1/2 simply by —61/2(0), because 
there is also an indirect effect of the changed global C(9) as well as the changed monopole 
that is being subtracted. These two are approximately constant if only a small area of the 
sky is affected, so analysis of 61/2 (&) is suitable for examination of the effect of small regions 
of the sky on Sx/2- However, these indirect effects must be taken into consideration while 
examining much larger regions. 



3. Directional contributions to 5*1/2 

3.1. Random map analysis 

We first examine plots of 5i/2 in random maps to obtain an idea of what to expect. 
We generate 80,000 synthetic Gaussian random statistically isotro pic maps, based on the 



underlying best-fit ACDM cosmological model. (ILarson et al.ll201ll ) 



A strong relationship between contributions to power in 61/2 and the regions of the 
quadrupole and octupole (QO regions for simplicity) is consistently seen in random maps. 
An example may be seen in the right panels of Figures [3] and H] with the full-sky ILC map. 

In order to quantify these contributions, we repeat our analysis of 61/2(0) in the regions 
of the sky containing only the quadrupole and octupole. We now compute 61/2 (O) in a 'sky' 
containing only the quadrupole and octupole of each map, and define the QO region to corre- 
spond to pixels Clgo where |5i/ 2qo (Oqo)| > x \ x is selected so that ^2 61/2 (Oqq)/ Yl ^1/2 (O) ~ 
0.9. We define the QO region of each individual synthetic map separately. The contribution 



2 For a masked sky, defining W(£li, 6) as the fraction of the ring for the respective C(Cli, 9) that is unmasked 
(W is if is a masked point) we have 

fi f 1/2 W(n t , 9)0(^,9)0(0)4008 6) 

17 J-i ELw(Ci v 9) 



Fig. 3. — The quantity Si/2(0) denned in Eq (jSJ) as a function of direction Cl in the cut-sky 
ILC map. The left panel shows contributions in the cut-sky map, where the masked pixels 
contribute nothing and are counted as zero. The right panel shows Si/ 2 (Cl) calculated from 
the full-sky ILC map. In both cases, the total value of Si/2(0) summed over the entire map 
is mathematically equal to 6*1/2, which is about 1000 (/iK) 4 for the cut-sky map, and 8000 
(/iK) 4 for the full-sky map. Note the strong resemblance of the full-sky S-\ /?.(&) and the 
power in the quadrupole and octupole (seen in e.g. Fig 4 in ICopi et al.l ( 120101 )). 




Fig. 4. — Si/2(&) is plotted in the cut-sky (left panel) and full-sky ILC map (right panel), 
with only the QO regions (defined in text) shown. Gray indicates both masked points and the 
non-QO region. In the full-sky map (right), correlation between the structure in £1/2(0) and 
the QO regions is apparent, as expected. In the cut-sky map (left), there is little correlation 
between the structure and the QO regions, a result at odds with typical statistically isotropic 
random maps. 
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of S1/2 from the QO region is thus defined as 

^ = E-M A ?°)- ( 7 ) 

i 

There is an extremely strong, almost linear relationship between S®j^ and total S1/2, 
plotted in Fig El In contrast, the total S1/2 and the contributions from non-QO regions of 
the sky are essentially uncorrelated, except at very high (> 10 5 (/1K) 4 ) values of S1/2, with 
which we are not concerned. Note that by definition, 6*1/2 = Sf, 2 + 5"°2 _Q °- 

This result makes intuitive sense, as the statistic Si/2{&) is based off of C(f2j, 6), which 
is higher when the central point Cl is in a high-temperature region such as the QO region. 
More qualitatively, S 1/2(^1) is a measure of correlation on the largest scales in the sky - which 
are, of course, the quadrupole and octupole. 



3.2. WMAP Map analysis 



We now examine WMAP's seven year maps in our analysis (IJarosik et al.ll201ll ). We 
use the full sky and cut-sky ILC (Internal Linear Combination) mapg^l The maps are first 
degraded to A^ si d c = 64 (corresponding to pixel scale ~ 1°); next, the KQ75 mask is applied, 
and finally all resulting pixels that are more than 10% masked (i.e. with mask value of less 
than 0.9) are excluded from the analysis. Finally, the dipole was removed. 

From inspection of the full-sky ILC S1/2 (fl) (Fig. HI right panel), there is an obvi- 
ous correlation between the magnitude of 5i/2(fi) contributions and the regions where the 
quadrupole and octupole have power, as expected from our analysis of random maps. Quan- 
tifying this result using our previous definition of the QO region, we find that 5^° = 8837 
(/1K) 4 , while S™- QO = -349 (/iK) 4 . The values observed in the full-sky map are very 
typical, compared to our random map results. 

Moving on to the cut-sky map, there is much less power on large angular scales, as 
expected; only 0.06% of random maps have a similarly low value of 5*1/2 . And as the empirical 
relationship predicts, 5^ is similarly low, at the 0.046% level. In comparison, the non-QO 



3 We have checked these results for robustness by comparing to a co-added foreground-cleaned Q-V-W 
band map, and the results there are very similar to that of the cut-sky ILC. 

4 Note that S^ 2 can be greater than the full S1/2, since the contributions from each individual S\/2 (Oqo) 
can be negative as well as positive, and their sum over the full sky is equal to the full Si/2- 
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contribution, Sy^'® is quite typical (15.6% confidence level.) Thus, the strong correlation 
that we observed in random maps continues to hold for the anomalous cut-sky map, as shown 
in Figure [21 What is unlikely about our sky is simply for S^® to exhibit such a low value in 
the first place. 

It is known that the quadrupol e in the CMB sky is lower than expected from ACDM 



results, though not anomalously so. ljKomatsu et al.l 120111 ) It makes sense to consider what 



happens to our results above if we use the ACDM theoretical power spectrum for Monte 
Carlo simulations as before, but substituting the low WMAP quadrupole (C 2 = 210 (fiK) 2 ), 
in order to test if our result is still valid. As expected, the significance of Si/ 2 relative to these 
new synthetic maps is now lower (0.70% compared to 0.06% before), and a similar effect is 
seen in our new statistics: Sy 2 has a significance of 0.58% (previously 0.046%), and 
is again typical - 27.4% C.L., compared to 15.6% before. Therefore, comparison to Gaussian 
random skies conditioned to have low C2 gives results that are somewhat less significant then 
before, which is expected since we already knew that part of the low-S'1/2 problem originated 
from low quadrupole. However, the strong correlation observed in random maps continues to 
hold. In addition, the various statistics are still low with < 1% probability, even compared 
to skies with the low quadrupole, and the lowness is once again due to lack of power from 
the quadrupole and octupole regions. 



4. Conclusion 

In this paper, we first introduced the fixed- vertex correlation function C(8,Cl), which 
has the potential to be a useful tool for studying the statistical isotropy of the correlation 
function. We then observe that the value of S1/2 in random skies is correlated with an 
expected anisotropy in the quadrupole and octupole regions of the sky. Our sky, known to 
have a low 5*1/2, also receives low contributions from its quadrupole and octupole regions. 
We have a bit of a chicken and egg problem here, as we cannot tell if the low Sy 2 caused a 
low S\/2i or vice versa, or if bo th were caused by an und erlying factor. It is worth considering 



that this result and the one of iPontzen fc Peirid (120101 ) may stem from the same root cause. 



It is also interesting that random Monte Carlo skies (which are generated from the 
assumption of statistical isotropy) and the reconstructed full-sky ILC (which assumes sta- 
tistical isotropy in the reconstruction process) have similar values of 5*1/2 and Sy®, while 
the cut-sky temperature map is anomalous. These findings may provide a hint for a future 
understanding of the origin of large- angle missing power in the CMB. 
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